DFTT 5/2001 
HUB-EP-00/05 
ITEP-TH-81/00 



Bound states and glueballs in 
three-dimensional Ising systems 

M. Caselle", M. Hasenbusch^, P. Provero'^'" and K. Zarembo'^'*^ 

" Dipartimento di Fisica Teorica dell'Universita di Torino and 
Istituto Nazionale di Fisica Nucleare, Sezione di Torino 
via P.Giuria 1, 1-10125 Torino, Italy 

^ Humboldt Universitat zu Berlin, Institut fiir Physik 
Invalidenstr. 110, D- 101 15 Berlin, Germany 

Dipartimento di Scienze e Tecnologie Avanzate 
Universita del Piemonte Orientale, Alessandria, Italy 

Department of Physics and Astronomy 
University of British Columbia, Vancouver, BC V6T IZl, Canada 

^ Institute of Theoretical and Experimental Physics, 
B. Cheremushkinskaya 25, 117259 Moscow, Russia 

E-mail addresses: 
caselleSto . inf n . it , hasenbusSphysik . hu-berlin . de 
proveroOto . inf n . it , zaremboOphysics . ubc . ca 

Abstract 

We study the spectrum of massive excitations in three-dimensional 
models belonging to the Ising universality class. By solving the Bethe- 
Salpeter equation for 3-D 0^ theory in the broken symmetry phase we 
show that recently found non-perturbative states can be interpreted 
as bound states of the fundamental excitation. We show that duality 
predicts an exact correspondence between the spectra of the Ising 
model in the broken symmetry phase and of the Z2 gauge theory 
in the confining phase. The interpretation of the glueball states of 
the gauge theory as bound states of the dual spin system allows us to 
explain the qualitative features of the glueball spectrum, in particular, 
its peculiar angular momentum dependence. 



1 Introduction 



The concept of universality is one of the fundamental principles in the theory 
of critical phenomena. According to the universality hypothesis, the long- 
range behavior of a statistical system near a phase transition is governed by 
fluctuations of order parameters and is essentially determined by symmetries, 
no matter how complicated the underlying microscopic dynamics is. For 
systems with Z2 symmetry universality leads to the description in terms of 
the one-component 0^ field theory (see e.g. Refs. [1, 2]). 

It is known that an order parameter in three-dimensional systems does not 
fluctuate strongly and the mean-field approximation (Landau theory) works 
reasonably well. The systematic methods that start from the mean-field 
approximation and perturbatively take into account fluctuative corrections 
give very accurate predictions for universal quantities such as critical expo- 
nents. However, a detailed numerical study of correlation functions in the 
scaling regime of various systems with Z2 symmetry revealed deviations from 
the mean-field approximation which are not too large numerically, but at the 
same time have no qualitative explanation within perturbation theory. These 
deviations arise in the hierarchy of length scales in the phase with broken Z2 
symmetry. 

The only macroscopic scale in the problem, of course, is the correlation 
length ^. A generic correlator decreases as e~^/^ at infinity, but there are 
some correlation functions that decrease more rapidly. Composite operators 
in the Ising model whose correlation functions decay faster than e~^/^ were 
constructed in [3, 4] and include, as typical examples, operators of non-zero 
angular momentum. What are the length scales associated with these opera- 
tors? The answer is simple in the Gaussian approximation, when correlators 
factorize: a correlation function of composite operators that factorizes on 
n irreducible correlators behaves as (e"''/^)". The characteristic decay rate 
of any correlation function is therefore an integer multiple of the inverse 
correlation length. 

The fact that fiuctuations are not exactly Gaussian should not change 
this conclusion. The ratio of the correlation length to any other length scale 
is an integer to any order in perturbation theory! This can be easily un- 
derstood from the field theory perspective. The behavior of a correlation 
function at infinity is determined by the singularities of its Fourier trans- 
form in the complex momentum plane. These singularities are associated 
with particles described by the quantum field. The only singularities that 
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arise in perturbative theory are the single-particle pole associated with 
the elementary quantum and the multiparticle cuts. Correlation functions 
of operators that couple to the single-particle state decay as e^^^, where /x 
is the inverse correlation length, or the particle mass. Correlation functions 
of operators that have zero overlap with the one-particle state (because of 
the angular momentum conservation, for instance) do not have poles in the 
complex momentum plane. Their only singularities are multiparticle branch 
cuts that start at p"^ = —n^ij,^. Consequently, such correlators decay as e~"''^^ 
with some integer n. 

The numerical simulations clearly show that the above picture does not 
hold in the phase with broken Z2 symmetry. The firmest evidence comes 
from the existence of a spin zero state with mass just below the two-particle 
threshold: M ~ l.Sfi, strongly supported by numerical results [3, 4, 5]. This 
and other non-perturbative states give visible contributions to the quanti- 
ties which are sensitive to distances of order of the correlation length. For 
instance, the spin-spin correlator in the Ising model near the critical point 
cannot be fitted by a simple pole-plus-cut structure suggested by perturba- 
tion theory and definitely contains an extra pole below the cut [3]. 

The simulations show that dimensionless ratios like M/ jj, become constant 
when approaching the critical temperature. Moreover they have the same 
value within statistical errors in two different realizations of the universality 
class: the Ising model itself and the lattice regularized (j)^ theory. These two 
facts allow us to be confident about the universal character of these effects. 
Specifically, we expect the same spectrum to appear in all realizations of 
the universality class, as long as one moves away from the critical point 
by perturbing with a Z2-symmetric operator towards the region where the 
symmetry is spontaneously broken. 

To resolve the aforementioned contradiction, we proposed to interpret the 
non-perturbative states in the broken symmetry phase of 0^ theory as bound 
states of elementary quanta of the scalar field [4]. We shall develop this idea 
further in this paper. 

Let us briefly recollect some facts about bound states. Two-particle 
bound states show up as poles in the Fourier transform of correlation func- 
tions below the branch cut associated with the two-particle threshold. In 
three dimensions, the threshold singularity is logarithmic, and, though the 
two-particle cut comes from one-loop Feynman diagrams and thus is accom- 
panied by a coupling constant, the smallness of the coupling is compensated 
by the large logarithm. In the vicinity of the threshold, perturbation the- 



2 



ory breaks down. Fortunately, the large logarithms can be resummed. This 
resummation substantially changes an analytic structure of correlation func- 
tions near the threshold and produces a pole that is invisible in perturbation 
theory. This pole is associated with a two-particle bound state. Its position is 
determined by the Bethe-Salpeter (BS) equation, which can be solved order 
by order in perturbation theory. 

In this work we analyze the BS equation for 3D (j)'^ theory in the broken 
symmetry phase. The analysis shows that a bound state of the elementary 
quanta indeed exists. We compute its binding energy at the next-to-leading 
order. While the value we obtain at the leading order is in very good agree- 
ment with the Monte Carlo data, the series appear to be badly behaved, so 
that a reliable numerical estimate of the subleading corrections to the bind- 
ing energy would require the knowledge of a large number of terms in the 
series. 

The 3D Z2 gauge theory is related to the Ising model by an exact duality 
transformation. Since the duality transformation maps the confining phase 
of the gauge theory into the broken symmetry phase of the spin model, all 
of our considerations apply to the glueball spectrum of the gauge theory as 
well. Indeed, we shall prove that duality implies the exact equality between 
the spectra of the spin model in the broken symmetry phase and the gauge 
theory in the confining phase for all values of the coupling. 

This correspondence suggests a novel interpretation of the glueballs of the 
gauge theory as bound states of the dual spin model, which in turn allows 
one to understand the qualitative features of the glueball spectrum, and in 
particular its peculiar angular momentum dependence. Bound states with 
different values of the angular momentum are easily seen to appear in the 
low-temperature expansion of the 3D Ising model. Simple group-theoretic 
arguments determine how many elementary quanta are necessary to construct 
a bound state of any given angular momentum and parity. Assuming that 
the binding energy is always much smaller than the mass of the elementary 
components, and that this description of the spectrum in terms of bound 
states is conserved when one moves from the deep low-temperature region to 
the scaling region, one ends up with a prediction for the qualitative features 
of the glueball spectrum in Z2 gauge theory, and in particular for the angular 
momentum dependence of the various masses. This prediction is in very good 
agreement with the numerical data for the glueball spectrum, thus lending 
support to the interpretation of glueballs as bound states of the dual spin 
model. 
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The rest of the paper is organized as follows: in Sec. 2 we will present the 
detailed computations for the bound state energy in the framework of the 
BS equation for 3D 0"^ field theory; in Sec. 3 we study the dependence of the 
binding energy on the magnetic field. Sec. 4 contains a proof of the exact 
correspondence between the spectrum of the 3D Ising model in the broken 
symmetry phase and the Z2 gauge theory in the confining phase. In Sec. 5 
we describe in detail the argument, already presented in [4], that allows us 
to explain the angular momentum dependence of the glueball spectrum of 
the Z2 gauge theory using the dual bound state picture. In Sec. 6 we discuss 
obtained results and outline possible future developments. 



2 Bound states in 3D 0^: BS equation 



We consider 3D 0^ theory in the broken symmetry phase. The Euchdean 
action is 



(2.1) 



The field 4> acquires an expectation value and the perturbative expansion is 
performed around the stable vacuum cf) — v. Extracting the fluctuating part, 



and setting 

we write the action as 



m 
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(2.2) 



(2.3) 



(2.4) 



The quartic coupling g has a dimension of mass, so the dimensionless expan- 
sion parameter is the ratio g/m. 

If the coupling is weak, the bound state lies very close to the threshold. 
We shall see that the binding energy is exponentially small in g/m. While 
the exponent can be easily found [4], computing the pre-exponential fac- 
tor requires sufficiently comphcated next-to-leading order calculations. The 
calculations are to certain extent similar to those of Ref. [6] where the BS 
equation was solved to the first few orders in perturbation theory for (1+1)- 
dimensional scalar field theory with the cosine interaction. 
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2.1 Leading order solution 




b 

Figure 1: a) Ladder diagrams, b) BS equation. 

The two-particle bound state manifests itself as a pole in the four-point 
function. The pole arises after summation of the ladder diagrams (Fig. la). 
Each ladder diagram develops a kinematic singularity when the momentum 
flowing through the diagram is near the two-particle threshold: ^ — 4m^, 
since the momentum integration in each loop then contains a region where 
two of the intermediate states simultaneously go on-shell and the propagators 
in the loop integral simultaneously have poles. This singularity produces 
large logarithms that compensate for powers of the coupling, so near the 
threshold all ladder diagrams are of the same order. The vertices in Fig. la 
correspond to the sum of the two-particle irreducible (2P1) diagrams, that 
is, the diagrams that cannot be split into disconnected pieces by cutting two 
internal lines whose total momentum is p. 2PI diagrams are analytic at the 
two-particle threshold in the s-channel. 

The position of the pole is determined by the BS equation, graphically 
represented in Fig. lb. The BS equation is the homogeneous part of the 
Dyson equation for ladder diagrams [7]: 



r(pi,P2) = 




T{q,p - q)G{q)K{q,p - q;Pi,P2)G{p -q), p = p^+p^. 



(2.5) 

Here, G is the propagator and K is the BS kernel, the aforementioned sum of 
2PI diagrams. The BS equation is a linear eigenvalue problem for the wave 
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function F and the bound state mass: 

= -p2- 

We take: 

M = (2 - A)m, 

where A is small. 



(2.6) 



(2.7) 



+ 



+ 



+ 




Figure 2: Tree-level BS kernel. 

To the leading order in g/m^ the BS kernel is given by the sum of three 
diagrams in Fig. 2. The diagram a contributes —g/2 to the kernel: 



9_ 
2' 



The diagram h gives 



3gm^ 



2 m2 - M2 



-| + 0(A). 



The contribution of the diagrams c and d depends on the momentum transfer 
in the t and u channels. However, the momentum transfer is very small when 

external momenta and both of the propagators in (2.5) are on shell: simple 
kinematic arguments show that the typical momentum in the t or u channel 
propagators is of order Am. To the first approximation, we can just set it to 
zero: 

3gm'^ 1 3g 

2 (pi - g)2 + m2 ^ T' 



3gm^ 



3^ 



2 (P2 - qy + m? 
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Collecting all three terms together, we get: 



K = 2g. (2.8) 

The kernel corresponds to a local four-point interaction in this approxima- 
tion. The effective quartic coupling appears to be negative. This happens 
because the exchange diagram c dominates over the pure four-point vertex a 
and the annihilation graph b, so the interaction is attractive. No matter how 
weak the interaction is, it will bind elementary quanta into bound states. 

Since the binding energy is small, the non-relativistic approximation 
should be valid at small g/m. The non-relativistic limit of the BS equa- 
tion with the positive local kernel is the two-body Schrodinger equation with 
attractive 5-function potential. This quantum-mechanical problem has a 
bound state at any value of the coupling with exponentially small binding 
energy [8]. An approach based on the non-relativistic approximation was 
elaborated in [4]. This approach is simple and physically transparent, but it 
only allows us to find the binding energy up to a numerical factor. Here we 
follow a more systematic route based on the full BS equation. 

The BS equation with the local kernel Eq. (2.8) corresponds to summing 
an infinite series of bubble diagrams, and can be written as 

_ f d^q 1 
^~^^J (2^ (?2 + ^2)[(p_^)2 + ^2]- (2.9) 

Evaluating the integral, 

(27r)3 (^2 + ^2-) _ ^^2 + ^2] 47r-v/p2 SttM 

^^^+0(AlnA), (2.10) 



we obtain: 



This equation predicts an exponentially small binding energy: 

A = const • exp i J , (2.12) 
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in accord with the result obtained in the non-relativistic approximation [4]*. 

The normahzation of the binding energy cannot be deduced from Eq. (2.11). 
We cannot keep the constant term in (2.10), because it has the same magni- 
tude as logarithmic terms in the next order of perturbation theory, which are 
proportional to gi^/m^ In A ~ gi/m. There are other sources of corrections of 
order g/m, as well. We need to compute all of them to fix the normalization 
of the binding energy. 

2.2 Corrections to propagator 



One-loop corrections to the propagator are given by three diagrams in 
Fig. 3. They shift the pole of the propagator to 




Figure 3: One-loop corrections to the propagator. 




3gm 



In 3 - ^A, 



(2.13) 



167r 



where by A we denote the linearly divergent integral 




(2.14) 



Expanding the third diagram in 



Fig. 3 near the mass shell, we get 



G-\k) = {p^ + fj,^) 1+ 3 



--ln3 



247rm 
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(2.15) 



and 





*The coupling A used in [4] differs from 5 by a factor of 24. 
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In evaluating the integral (2.10), the corrections to the propagator will result 
in the replacement of the bare mass m by the renormalized (physical) mass 
fi and in the multiplication of the whole integrand by the wave function 
renormalization factor. 

It is natural to express everything in terms of the physical mass gap 
(2.13). In particular, we set: 



M = (2 - A)/i. 



(2.17) 



This definition affects the BS kernel, since the contribution of the diagram b 
in Fig. 2 depends on how exactly we define M: 



3gm^ 



1 



3^ 



m 



9 
2 



,2 



9 



In 3 



2 - M2 

2g'A 



1 - 



3gm 



167r 

+ 0(A) + 0(^3). 



In 3 + 



(2.18) 



Snm 3m2 

This differs from the result based on our previous definition by an amount 

^2 In 3 2g^A 



87rm 3m2 

2.3 One-loop corrections to the kernel 



(2.19) 




1 





7 8 9 10 11 12 13 

Figure 4: One-loop two-particle-irreducible diagrams. 

The one-loop corrections to the kernel of the BS equations are given by 
thirteen diagrams in Fig. 4. Since we are interested only in the logarithmically 
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enhanced terms in the BS equation, all these diagrams can be evaluated on 
shell and at zero momentum transfer in the t and u channels. 
The linearly divergent diagrams 1-3 sum up to zero: 

5Xi_3 = 0. (2.20) 

The diagrams 4-7 cancel the divergence in the diagram b, the second term in 
Eq. (2.19): 

«,.. = |^. (2.21) 

All linear divergences eventually have canceled. It is satisfying to see that 
the net one-loop contribution to the kernel is finite, as it should be. 
The remaining diagrams give: 

5^8-13 = (9-6 + 1 + 12-4 + 2)-^ = ?-^. (2.22) 

Collecting together the contributions (2.19)-(2.22) we get the one-loop cor- 
rection to the BS kernel: 

6K = (7 -InS)-^. (2.23) 



2.4 Momentum dependence 

Finally, we consider the most subtle corrections that come from the momen- 
tum dependence of the kernel and of the wave function. Those were neglected 
in the leading order calculation because the logarithmic enhancement of the 
bubble diagram near threshold comes from the region of integration in (2.5) 
where all momenta arc almost on shell. This allowed us to forget about 
the momentum dependence of the wave function and to set it to one. In 
fact, we can always require that the wave function is equal to the one on 
shell. This can be regarded as a normahzation condition. But to compute 
non-logarithmic terms we need to know the complete off-shell wave function. 

Once the on-shell condition on the external momenta pi,P2 is relaxed, 
we can no longer neglect the momentum transfer in the diagrams c and d in 



10 



Fig. 2. The BS equation then gives: 



r(pi,P2) 



9 



Sq 1 



?>grn? 



{q + m2 



87rm 



ln( ^ 



4m^ + Pi + pI 
1 



4777.2 2 \ 777/ 

2.5 Next-to-leading order BS equation 

With all corrections taken into account, the BS equation reads: 



(2.24) 



1 = 



X 



X 



3- 

d^q 



9 



- InS , 

4 / 247r/i 



(27r)3 

-Pl)^ + 777^ 



4/x2+p2_^(p-g)2 



g + {7-]n3) 



(2.25) 

where external momenta pi and p2 are on shell: p^ 2 = ~A*^- We need to 
retain terms of order g/n and {g'^ / jji^) In A on the right hand side. For this 
reason, 0{g/m) corrections to the wave function can be neglected, since they 
vanish on shell and therefore do not lead to the logarithmic enhancement. 
The loop integrals encountered in (2.25) are listed in Appendix and in (2.10). 
Using these results, we get after a little of algebra: 



log A 



■^ + 21og2- Jlog3 + 
9 2 



(2.26) 



that is''' 



A = — = exp 



9 



(2.27) 



'''Note that the pre-exponential factor quoted in the published version of Ref. [4] is not 
correct. 
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2.6 Comparison with numerical results 

To compare Eq. (2.26) with numerical results one has to know a value for 
the dimensionless couphng u = g/ji. The safest thing to do is to relate u to 
the renormalized coupling ur defined e.g. as in Ref. [18], whose critical value 
is precisely known from Monte Carlo simulations [13] 



Ur = 14.3(1) 

in terms of ur we have, using results from [3, 18] 



u 



Ur 



Ur 



/ 15 3 log 3 
Vl287r ^ 327r 



so that 



Stt / 15 log 3 ^, ^ , ^, , 

logA = + +^ + 21og2 +OM 

Ur vis 4 



Therefore the leading-order result is 



logA = -— = -1.76(2) 

Ur 



corresponding to 



(2.28) 



(2.29) 



(2.30) 



(2.31) 



A = 0.172(3). 



(2.32) 



This result is compatible with the assumption of weak binding energy that 
underlies the Bethe-Salpeter formalism, and is in a very good agreement with 
the numerical result. 



Amc = 0.17(3). 



(2.33) 



However, the parameter of expansion in perturbation theory is comparatively 
large, and the NLO correction is not small: the 0(1) contribution to log A is 
~ 2.6, that is actually greater than the LO contribution. Therefore a rehable 
estimate of the subleading corrections to the binding energy seems to require 
the computation of a large number of terms in perturbative series, so as to 
make a resummation possible, which is typical for application of "id (J)'^ theory 
to critical phenomena. 
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3 Dependence on the magnetic field 



It is possible to study the dependence of the bound state masses on the 
magnetic field. Because the magnetic field couples to the order parameter, 
we know how to introduce it in the effective description in terms of (p^ theory. 
We can take the magnetic field into account by adding a linear term to the 
action: 



(Fx 



2\2 



h(t) 



(3.34) 



The magnetic field shifts the minimum of the potential, which now is deter- 
mined by the equation 



Expanding the order parameter around 0o, 
we get for the action: 



-I 



(fx 



+ So, 



where the mass and the trilinear coupling are 

9_ 
6 



(3.35) 



(3.36) 



(3.37) 



(3.38) 



93^ g 00- (3.39) 

The sum of the diagrams in Fig. 2 gives the tree-level BS kernel: 

= 9l .91 1 g _^ 7(l>l + v^ 

" ml ^ 2 ml- Ami 2 ^602- 27;^ " ^"'■^''^ 

Repeating the same steps as we used to calculate the binding energy at zero 
magnetic field, wc get 

A = const exp (^-^^^^^ , (3.41) 
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where A is the dimensionless binding energy: — (2 — A)mh- 
It is useful to introduce the dimensionless variable 




(3.42) 



which satisfies the equation 




(3.43) 



The binding energy is expressed in terms of ^ and m, the mass at zero 
magnetic field, as 



It is easy to sec that A is a decreasing function of h. Thus the magnetic 
field loosens the binding of composite states and shifts their masses closer to 
the threshold. 

4 Duality and the spectrum of Ising systems 

The 3D Ising model and Z2 gauge theory are related by an exact duality 
transformation. The broken symmetry phase of the spin model is mapped 
into the confining phase of the gauge theory. The purpose of this section is to 
show that, in this phase, the duality relationship implies an exact coincidence 
of the spectra of the two theories. The proof relies on the existence of a 
non-zero interface tension and therefore applies exclusively to the broken 
symmetry phase of the spin model. 

There are several books and reviews which discuss duality in spin systems 
and in particular in Ising- type models (see e.g. Ref. [9]). However the duality 
transformation is usually treated in the thermodynamic limit only, while 
to study the spectrum of the transfer matrix it is necessary to extend the 
analysis to finite lattices. This is possible as long as one carefully takes into 
account all possible boundary conditions. 

The novel feature of our approach, which greatly simplifies the whole 
analysis, is the use of the Transfer Matrix (TM) formahsm. Let us first 
recall the definitions of the lattice models we are concerned with, and the 
notion of duality in the thermodynamic limit. 




(3.44) 
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The spin Ising model 

The Ising model is defined by the action 

^spin f^spin ^ ^ ^n^m i (4.45) 

(n,m) 

where the field variable s„ takes the values —1 and +1; n = (no, ^i, 712) 
labels the sites of a simple cubic lattice of size Lq, Li and L2 in the 
three directions. The notation {n, m) in Eq. (4.45) indicates that the 
sum is taken over pairs of nearest neighbor sites only. 

The gauge Ising model 

The building blocks of the Z2 gauge model are the hnk variables gn-,j, £ 
{ — 1, 1}, which play the role of gauge fields. Denoting by ji the direction 
of the link, the action is 

Sgauge = "/^ ^ Qn^tiv (4.46) 
n,ix<v 

where Qn-^fiv are the plaquette variables, defined by 

9n;^v 9n;^ 9n+fi;u 9n+v,fi 9n;v ■ (4.47) 

This action is invariant under local Z2 gauge transformations defined 
as follows: one chooses arbitrarily a subset of the sites of the lattice and 
changes signs of all variables defined on the links which end on these 
sites (if two neighboring sites belong to the chosen subset the hnk that 
joins them is changed twice, that is not changed). It is immediately 
clear that the plaquette values, like any other product of links along a 
closed path, are invariant under this gauge transformation. For more 
details about this model, see e.g. Ref. [10]. 

Duality 

There is an exact duality transformation which relates the Ising model 
and the Z2 gauge model. This transformation is known as Kramers- 
Wannier duality. It relates the partition functions of the two models 
evaluated at two different values of the coupling constants: 



Sgauge ( /3 ) CX! Zigpi^i ( /3 ) 



~P = -^log[tanh(/3)] , (4.48) 
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where (3 will be denoted as the "dual coupling" in what follows. 

It is easy to see that low values of (3 are mapped into high values of /9 

and vice versa. Thus the confining region of the gauge theory is mapped 
into the broken symmetry phase of the spin model. In particular the 
end points of these two phases, the deconfinement transition and the 
magnetization transition, are mapped into each other. 

An important feature of the dual transformation on a lattice of finite size 
is that it does not conserve the boundary conditions (BC). In the thermo- 
dynamic limit this fact becomes irrelevant, but on lattices of finite extent it 
cannot be neglected. In particular the Z2 gauge model with periodic BC in 
all directions is mapped by duality into the Ising spin model with fluctuat- 
ing BC, so that the partition function is given by the sum of the partition 
functions with all possible choices of periodic (p) or antiperiodic (a) BC: 

^ gauge,ppp(,P^ — (-^spm,ppp(/3) ~l~ Zgpiji^ppal^P^ -\- ... -\- •^spm,aao(^) )) (4.49) 

where V = LqLiL2 is the volume of the lattice, and c is a constant which 
can be easily evaluated, but is irrelevant for our purposes. 

This result is discussed in full generality (for a generic lattice geometry 
and symmetry group) in [11]. In the particular case of the Ising model on 
a cubic lattice it can be easily obtained by a direct implementation of the 
duality transformation. 

Let us now consider the duality transformation in the framework of the 
transfer matrix approach. The direction will be our "time" . If we choose 
periodic BC in the direction without specifying the BC in the 1 and 2 
directions we obtain: 

^spin,pxy tr(Txy) (4.50) 

where T^.^ denotes the transfer matrix of the model in which x,y & {a, p} 
BC are chosen in the 1 and 2 directions respectively. 

The antiperiodic BC in the direction can be obtained by acting with a 
spin-fiip operator P, which changes the sign of all spins in a given time slice. 
Thus we may write 

Since the operators T^y and P commute, they have a common set of 
eigenf unctions. Let us denote the eigenvalues of T^y by X^y^i and those of 
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P hy Pi. The possible values of Pi are 1 and —1. States that are symmetric 
in the magnetization have Pi — 1 and those that are antisymmetric have 

Pi = -1- 

Thus we can write 

Zspin,pxy ~l~ Zspin,axy — tr(l + P) (T^jj^) ° = 2 ^ ^ ^xy,ii (^■^^) 

i,sym 

where means that the sum is restricted to the states that are sym- 

metric in the magnetization (See Ref. [12] for further details on this type of 
construction). 

Using this result we can write the duality relation in terms of the transfer 
matrix eigenvalues: 



Lq „LqLiL2 



iySym i,sym i,sym i,sym 



, (4.53) 



where 7^- are the eigenvalues of the transfer matrix of the gauge system. 

It is instructive to show explicitly that the sums on the two sides of 
Eq. (4.53) have the same number of terms. On the spin side (right hand side 
of Eq. (4.53)) we have 2^^^^~^ terms for each sector. In fact, the transfer 
matrix is a 2^^^^ x 2^^^^ matrix, but only half of the eigenvalues fulfil the 
symmetry requirement. Thus taking together all the four sectors we end up 
with a sum of 2^1-^2+1 terms. On the gauge side we have 2L1L2 link variables, 
but we may fix the gauge. The maximum number of gauge variables that 
we may fix, in order to preserve the periodic BC, is exactly L1L2 — 1. This 
is easy to see: we may fix, say, all the links in the 1 direction except those 
in the last row, where we may fix all the links in the 2 direction, except the 
last one. Thus we end up with only L1L2 + 1 degrees of freedom left. Hence 
the corresponding transfer matrix has 2^^^^~^^ eigenvalues, that is exactly the 
same number as for the spin model. 

Since for finite values of Li and L2 the two sums contain a finite number 
of terms and since Eq. (4.53) holds for any integer Lq, each 7^ has to have 
an exact (up to an overall factor c^'^^'-^) counterpart in the left hand side of 
the equation. The overall factor cancels when one takes the ratios of all the 
eigenvalues to the lowest one to obtain the physical spectrum, and so the 
latter coincides in the two models. 

However our goal is to compare the spectra of the gauge theory and 
the spin model both with periodic BC in all directions, while Eq. (4.53) 
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involves a sum over many different choices of the spin model BC. Therefore 
we must find a relationship between the App j eigenvalues and those belonging 
to the other three sectors. This is easily done by noticing that in the low 
temperature phase of the Ising model with antiperiodic boundary conditions 
at least one interface has to be created in the system. Therefore at leading 
order (neglecting the cases in which more than one interface appear): 

Zppa ~ exp(-(T LqLi) Zppp, (4.54) 

where a is the interface tension. 

For the transfer matrix eigenvalues this gives 

\p,max ~ exp(crLi) Xpa,max (4.55) 

where, with the notation Xpp^max we denote the low lying states of the spec- 
trum-f . 

So we can conclude that for sufficiently large Li and L2, the low lying 
spectrum of the Z2 gauge theory in the confining phase with periodic bound- 
ary conditions coincides with the symmetric sector of the Ising spin model 
spectrum in the low temperature phase with periodic boundary conditions 
only. (This last argument was already presented in [13]). 

5 The dual bound state picture of glueballs 
in Z2 gauge theory 

In this section we show that the interpretation of the spectrum of three- 
dimensional models in the universality class of the 3D Ising model as bound 
states of the elementary quanta provides an analytical tool to explain the 
qualitative features of the glucball spectrum of Z2 gauge theory, and espe- 
cially its angular momentum dependence. 

Consider first the spin Ising model in the low-temperature region. Here 
the existence of bound states can be immediately inferred from the diagram- 
matics of the low-temperature expansion. This was shown in [14] for the 
4D case, where however all bound states are expected to disappear in the 
continuum limit due to triviality. 

■''Notice that this last argument does not hold in the 2D case, where the anti-periodic 
eigenvalues are only suppressed by a constant, independent of the system size. 
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The mechanism is very simple, and is best explained in the transfer ma- 
trix formalism. As discussed in [15, 16], to the leading order in the low- 
temperature expansion of the transfer matrix, all time slices are forced to 
have the same configuration of spins, and the the eigenvectors are given by 
all the possible spin configurations in a single time shce. The eigenvalues 
of the Hamiltonian (that is minus the logarithm of the transfer matrix) are 
given by the energy of each configuration. 

The ground state then corresponds to the configuration in which all the 
spins on any time slice are parallel; the first excited state corresponds to one 
flipped spin, which requires four bonds to be frustrated. To proceed, one has 
to flip two spins: if these are chosen in non-neighboring sites, the energy is 
just twice the one of the first excitation, since eight bonds will be frustrated. 
However one can flip two neighboring spins at the cost of frustrating six 
bonds only: the corresponding state is a bound state of the fundamental 
excitation, with mass just below the two-particle threshold. 

One can also construct states of given angular momentum and parity 
by choosing linear combinations of time slice configurations with nontrivial 
transformation properties under spatial rotations and parity refiections. This 
can be done systematically by using standard group theory results, namely 
the theory of representations of the dihedral group D^, which is the relevant 
group on a square lattice. This analysis is performed in Ref . [17] , to which 
we refer for details. 

The important point is that any given angular momentum requires a cer- 
tain minimum number Uc of spin flips. This means that a bound state of such 
angular momentum will be composed of at least ric elementary quanta. In 
Tab. 5 we have reported the values of ric for the various values of the angular 
momentum. For example. Fig. 5 shows the simplest lattice operator corre- 
sponding to angular momentum J = 2 and parity P = — 1, which requires 
ric = 3 elementary quanta, while Fig. 6 shows how states with quantum 
numbers J — 0,2 and P — +,— can be constructed with four elementary 
quanta. If we assume that the binding energy is always much smaller that 
the common mass of the elementary constituents of the bound states, we end 
up with a prediction for the angular momentum dependence of the glueball 
spectrum of the Z2 gauge theory in the strong coupling regime. 

It is not at all guaranteed, of course, that the spectrum will retain these 
same features when going from the strong coupling regime to the scaling 
region: however the data from Monte Carlo simulations of the gauge theory 
show that this is actually the case. The third column of Tab. 5 contains the 
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Figure 5: Three elementary quanta are required to construct a state with quantum 
numbers 2~. Dots represent flipped spins, sohd hnes are frustrated hnks, dashed 
hues satisfied hnks. 

Monte Carlo results for the glueball masses (from Ref. [17]), normalized to 
the mass of the lowest glueball. This shows that our dual bound state picture 
indeed explains the pecuhar qualitative features of the glueball spectrum. 

In particular, some characteristic approximate degeneracies in the spec- 
trum find a natural explanation in this picture: for example the states O"*"'** 
and 2^ are nearly degenerate because they are both bound states of nc = 3 
elementary constituents. Bound states of = 4 constituents give rise to 
the 0~ and 2^'* states, again explaining their near degeneracy, and the same 
applies to 0~'* and 1^ when interpreted as rzc = 5 bound states^. We do not 
expect these degeneracies to be exact, since there is no reason to expect the 
binding energies to be exactly the same for bound states of different angular 
momentum. 

Once the minimum number Uc of constituents for a given value of the 
angular momentum has been reached, it is easy to see that states with the 
same angular momentum can be constructed out of any number n > ric oi 
constituents. This suggests that the approximate degeneracies just discussed 
appears between pairs of states only because, in general, only two states can 
be detected with numerical methods in each channel. It is very likely that 

§Thc 2+ state could be realized with ric — 2, but a general theorem (see Ref. [17]) 
forces all states with J 7^ to be degenerate in parity. 
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Figure 6: Bound states of nc = 4 elementary quanta. 



the degeneracies actually group together larger and larger sets of states as 
the number of constituents is increased. For example, we might conjecture 
that the states 0~'* and 1^ {ric = 5) are nearly degenerate with the (as yet 
undetected) states 0+'****, 2±'**. This de generacy pattern is a prediction of 
the dual bound state picture of glueballs. 

6 Summary and further developments 

The broken symmetry phase of 3D models with Z2 symmetry shows a rich 
spectrum of massive excitations, contrary to naive expectations based on the 
perturbative field-theoretical treatment. The non-perturbative states in the 
spectrum can be interpreted as bound composites of the elementary particle 
excitation. When applied to 3D Z2 gauge theory, this approach gives a nice 
explanation of the angular momentum dependence of the glueball spectrum. 
In particular, an observed rather peculiar degeneracy pattern naturally arises 
in the bound-state interpretation. 

It would be interesting to consider the contribution of the bound states 
to universal amplitude ratios. The obvious way to do that is to use the lad- 
der approximation for correlation functions, which incorporates all diagrams 
that blow up near the two-particle threshold. It would be also interesting 
to understand if there are bound states in the systems with > 1 order 
parameters. The physics is completely different in this case because the 
dominating forces are long-range due to the Goldstone bosons of the broken 
0{N) symmetry. 
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Table 1: The bound states that can be constructed with tIc elementary quanta 
are classified in terms of their angular momentum and parity. The last col- 
umn shows the mass of the corresponding glueball state in Z2 gauge theory. 



Tie 




m(J^)/m(0+) 


1 


0+ 


1 


2 


0+* 


1.88(2) 


3 


g+** 


2.59(4) 




2± 


2.59(4) 


4 


2±* 


3.23(7) 




0" 


3.24(16) 


5 


Q-* 


4.48(20) 




1± 


4.12(17) 
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Appendix A 



d^q 1 

(27r)3 [4/^2 + q-^ + {p- 5)2] (^2 + ^2)[(p _ + ^2] - ^277 \9A 

(A.l) 



d^q 1 

{2t:Y [q^ + m^)[{p — qY + rn^][{q — PiY + rn?] 167r/i^ 

(A.2) 



/ 



(27r)3 



[Ajj? -\- q^ -\- {p — qY\{q^ + rn?)[{p — q^ + ni?][(q — pi)^ -|- m^] 
^ • ' — - (A.3) 



327r/i5 



i"(4V- 
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All integrals are calculated on shell: = —fJ''^, and near the two-particle 
threshold: = (p^ +P2f = -M^ = (2 - A)V^- 



References 

[1] C. Itzykson and J. M. Drouffe, Statistical Field Theory, (University 
Press, Cambridge, 1989). 

[2] J. Zinn- Justin, Quantum Field Theory And Critical Phenomena, 
(Clarendon, Oxford, 1989). 

[3] M. Caselle, M. Hasenbusch and P. Provero, Nucl. Phys. B556 (1999) 
575 [hep-lat/9903011]. 

[4] M. Cascllc, M. Hasenbusch, P. Provero and K. Zarembo, Phys. Rev. 
D62 (2000) 017901 [hep-th/0001181]. 

[5] D. Lee, N. Salwen and M. Windoloski, hep-lat/0010039. 

[6] R. F. Dashcn, B. Hasslacher and A. Ncveu, Phys. Rev. Dll (1975) 3424. 

[7] V. B. Berestetsky, E. M. Lifshits and L. P. Pitaevsky, "Quantum Elec- 
trodynamics", Pergamon, Oxford 1982. 

[8] C. Thorn, Phys. Rev. D19 (1979) 639. 

[9] R. Savit, Rev. Mod. Phys. 52 (1980) 453. 

[10] J. Drouffe and J. Zuber, Phys. Rept. 102 (1983) 1. 

[11] C.Grubcr, A.Hintcrmann and D.Merlini, Group Analyisis of Classical 
Lattice Systems, ( Springer- Verlag, Berlin-Heidelberg- New York, 1977). 

[12] M. Caselle, R. Fiore, F. Ghozzi, M. Hasenbusch, K. Pinn and S. Vinti, 
Nucl. Phys. B432 (1994) 590 [hep-lat/9407002]. 

[13] M. Caselle and M. Hasenbusch, Nucl. Phys. B470 (1996) 435 [hep- 
lat/9511015]. 

[14] M. Liischer and P. Weisz, Nucl. Phys. B295 (1988) 65. 

[15] M. E. Fisher and W. J. Camp, Phys. Rev. Lett. 26 (1971) 565. 

23 



[16] W. J. Camp, Phys. Rev. B7 (1973) 3187. 

[17] V. Agostini, G. Carlino, M. Caselle and M. Hasenbusch, Nucl. Phys. 
B484 (1997) 331 [hep-lat/9607029]. 

[18] G. Miinster and J. Heitger, Nucl. Phys. B424 (1994) 582 [hep- 
lat/9402017]. 



24 



